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Abstract 

We compute one-point function of the glueball loop operator in the maximally 
confining phase of supersymmetric gauge theory using instanton calculus. In the 
maximally confining phase the residual symmetry is the diagonal U(l) subgroup 
and the localization formula implies that the chiral correlation functions are the 
sum of the contributions from each fixed point labeled by the Young diagram. 
The summation can be performed exactly by operator formalism of free fermions, 
which also featured in the equivariant Gromov-Witten theory of P 1 . By taking the 
Laplace transformation of the glueball loop operator, we find an exact agreement 
with the previous results for the generating function (resolvent) of the glueball 
one-point functions. 
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1 Introduction 



Symmetry helps us to understand physics, in particular, the non-perturbative dynamics. 
Guided by holomorphy coming from M = 2 supersymmetry and the asymptotic behavior 
in semiclassical region, Seiberg and Witten wrote down the effective prepotential of M = 2 
super symmetric gauge theory pQ. However, it is still desirable to derive their result following 
the standard process of quantum field theory - path integral. After a lot of works on direct 
integration over the moduli space of instantons, Nekrasov [2] applied the localization theorem 
for toric actions to compute the equivariant integral over the instanton moduli space. Subse- 
quently, it was noted by [3] that a similar approach also applies to the instanton calculus of 
N = 1 supersymmetric gauge theory where J\f = 2 theory is perturbed by a superpotential 
Here the integration is localized to the summation over the (isolated) fixed points 
labeled by the Young diagrams, which are in some sense regarded as saddle points. One might 
expect the instanton calculus on the saddle points has more applications than kinematical 
constraints. 

The computation of the holomorphic quantities in M = 1 theory such as effective super- 
potential is related to matrix model in [3]. This proposal is derived by relating them to the 
B-model topological string amplitudes on certain non-compact Calabi-Yau manifold. The re- 
lation is a mirror to a gauge/string duality due to the geometric transition on the A- model side 
[5] . Thus the reproduction of the holomorphic quantities in M = 1 theory from the microscopic 
instanton calculus can be regarded as a check of the gauge/string correspondence involving 
the matrix model. For a recent progress in microscopic approach to M = 1 supersymmetric 
gauge theory, see [6j [7J E] • 

Though the localization formulas are given explicitly for both the partition function and the 
chiral correlation functions, in general it is not straightforward to perform the summation over 
the Young diagrams exactly. In the maximally confining phase of U(N) supersymmetric gauge 
theorjQ, where SU(N) is confined and the residual symmetry at low energy is the diagonal 
U(l) C U(N), the task of summation gets considerably tractable. In [9], we adopted the 
standard correspondence between Young diagrams and neutral states in the fermion Fock space 
and calculated explicitly the chiral one-point functions (Tnp J ) in the maximally confining 
phase. We found that the results can be summarized compactly in terms of the loop operator: 

<Tre^> = I U/Q ^f^ ) - VV*") , (ft - 0) , (1.1) 

with I n (x) being the modified Bessel function. The vacuum expectation value of the loop 
operator has two parameters q = A 2N and h. The series expansion in q is the instanton 

1 We reduce the computation to the instanton calculus of U(l) gauge theory. However, the existence of 
underlying U(N) theory should be assumed, since we will consider the expansion in the second Chern number 
whose meaning is lost in a genuine U(l) theory. 
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expansion of gauge theory, while that in H is interpreted as the genus expansion of the corre- 
sponding string theory. In the gauge/string theory correspondence, the partition function and 
the correlation functions allow the double perturbative expansion in the instanton number of 
gauge theory and the genus of the string theory. Usually we can sum up the expansion only 
in one of the two expansion parameters, but have to compute the functions order by order in 
the other parameter. It is very amusing that the above result (11.11) gives an example where 
we can perform the summation of the double expansion completely in a closed form. 

The chiral one-point functions (Tr ip J ) in Af = 1 theory do not detect the superpotential 
which is used to deform the original M = 2 theory [3J, [9] . In this paper we would 
like to proceed to the computation of the chiral glueball one-point functions (Tr X a X a ip J ), 
which depend on the details of W($>). The chiral glueball one-point functions were obtained 
previously by several methods. For example, one can relate the holomorphic (F-term) quan- 
tities in super symmetric gauge theory to amplitudes of topological string theory |IQj and then 
compute with matrix model using open/closed string duality [1]. We can also relate the cor- 
relation functions to the generalized Konishi anomaly and perform a purely field theoretical 
computation [TT]. The resulting glueball resolvent is given by 



with an unknown polynomial f(z) which is determined so that the asymptotic behavior of 
R{z) in z — > oo is (Tr \ a \ a ) ■ z^ 1 + 0(z~ 2 ). In the maximally confining phase, the polynomial 
f(z) is further tuned so that the square root is factorized into 



In matrix model the maximally confining phase is described by a one-cut solution, while on 
the moduli space of Seiberg-Witten theory it corresponds to the degenerating loci where the 
maximal number of mutually non-intersecting cycles of the Seiberg-Witten curve collapse. In 
this paper we would like to add one more method to derive the glueball resolvent R(z). We 
find our instanton calculus gives an exact agreement to the result (11. 2p with the factorization 
(11.31) . The reproduction of previously known results should serve as a non-trivial consistency 
check of the instanton calculus. Besides, along the way of computation we also find an explicit 
expression of the polynomial f(z) that causes the complete factorization ( 11. 3ft . It is interesting 
that our £7(1) instanton calculus automatically gives a closed form for the glueball resolvent 
in the maximally confining phase. 

The plan of the present paper is as follows. We will shortly review some necessary ingredi- 
ents of the instanton calculus in section 2. Before proceeding to our computation of the chiral 
glueball one-point functions in section 4, we will first compute two-point functions of the loop 
operators as a preparation in section 3. Finally in section 5 we comment on the higher genus 
correction to our computation. 
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(1.2) 




(1.3) 
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2 Review of instanton calculus 



Chiral operators O in super symmetric field theories are, by definition, those annihilated by 
the fermionic charges Q h of one chirality. Two chiral operators are defined to be equivalent, if 
the difference is Q A -exact. The set of chiral operators is closed under the multiplication and 
form a ring, which we call chiral ring. From the supersymmetry algebra in four dimensions, 
{Qa, Qa\ = a aaP^ we can see that correlation functions of chiral operators are "topological" 
in the sense that they are independent of the positions of operators. Especially, topological 
one-point functions characterize the phase structure of vacua. In the four dimensional M = 1 
supersymmetric gauge theory with a single adjoint matter $, we have a vector multiplet 
(Afj,, X a ) and a chiral multiplet $ = (ip,ip a ), which are all in the adjoint representation. One 
can show that the generators of chiral ring are of the form Tr<p J , Tr X a <p J and Tr X a X a (p J 
The correlation functions of the chiral operators O are defined by 

{0) M=l := -J— III 0)exp(-% =1 ) , (2.1) 



VZ M 

=1 JM kjc 2 ) 

where the action for Af = 1 supersymmetric gauge theory SV=i is that of Af = 2 theory SV=2 
perturbed by a superpotential W(&): 

SV=i := 5^- =2 + J dx A d6 2 W($) . (2.2) 

The correlator {0)j^=i is normalized by the volume V of the non-commutative C 2 and the 
partition function Ztf=\ := J M exp(—Sjy=i). The integral is over the moduli space M. of U(N) 
instantons described, for example, by the ADHM construction. The moduli space has the 
decomposition M. := U k A4pf ) k, that is, the correlation function is a sum over the contributions 
from each moduli space A4n^ with a specific instanton number k. The correlation function 
(0)j\f = i for O = Ti(p J or Tr X a X a (p J is of our prime interest in this paper. We shall omit the 
subscripts Af = 1 or A/" = 2 of the correlators hereafter as long as it is clear from the context. 

The strategy to calculate these correlation functions is to first replace all the operators by 
their equivariant extensions [3]: 

Tr</^a (2j2) ATrj^ 7 , (2.3) 

Tr A ° A ^ J " - (7TWTI) aM A Tr pM ' (2A) 

d 2 6 W($) i-> a ( 2,o) A Tr W{F) , (2.5) 
with a being the equivariantly closed forms 



«(o,o) 

«(2,0) 
a (0,2) 
"(2,2) 



1 , (2.6) 

dz 1 A dz 2 + ihz l z 2 , (2.7) 

dz 1 A dz 2 - ihz 1 ^ 2 , (2.8) 

"(2,0) A (2(0,2) • (2.9) 
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The curvature T of the universal bundle over C 2 x M. N,k is expanded according to the direct 
product structure of the base space as follow^: 

F =: F + ty + ip 

= F^dx»dx v + {\ m dz m + tp m dz m } + {{d M U^){d M U) - U ] C^U} . (2.10) 

If we untwist the theorjH, the components F^ v and A m are identified with the field strength 
and the gaugino, which gives a vector multiplet in M = 1 theory. The remaining pair ijj m and 
(p gives a chiral multiplet. Then due to the localization theorem, our computation reduces to 
picking up the value at the fixed points. Since the fixed points are classified by the Young 
diagrams, the result is given in terms of summation over the Young diagrams. Especially in 
the maximally confining phase, the correlation functions (12. ip simplifies into 

where the partition function is given by 



k=0 J \Y\=k 



z ~T.l En (ft(D)) r ( 2 - 12 ) 



We have introduced the parameter q of the instanton expansion, h is the equivariant parameter 
associated with the U(l) action T h : (zi,z 2 ) — > {e %h Z\,e~ lh z<i) on C 2 . Note that the partition 
function Z depends on these parameters only through the combination q/h 2 . We denote by 
\Y\ the total number of boxes of a Young diagram Y. The hook length at a box □ G Y is 
denoted by h(ri) and the weight nney(^( D ))~ 2 * s called the Plancherel measure. At each fixed 
point Y we can estimate the chiral operator O to obtain Oy. Thus Oy is a function on the 
space of Young diagrams and (O) is nothing but the integration of Oy with respect to the 
Plancherel measuraj. 

For the loop operator O = Tre tv , which is a generating function of Tr (p J , the function Oy 
is given by 

Ch y (a) = e ta (l + sh 2 (a) ^ e a(c(D)) ) , (2.13) 

OEY 

where a = (Try?), a = th and sh(a) := e Q / 2 — e~ a ^ 2 [9]. For a box □ at the m-th row and the 
n-th column of the Young diagram, we define the content by c(q) := n — m. In [3] the glueball 



2 The matrix U appears in the ADHM construction as the zero modes of the Dirac operator. 

3 Note that we consider the case C 2 ~ K 4 in this paper. In general the topological twist of Af = 1 theory is 
possible on the Kahler manifold. 

4 The Plancherel measure can be regarded as a discretization of the Vandermonde measure. This interpre- 
tation suggests a natural connection to the matrix model. 
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operator O = Tr X a X a (p J is related to the connected two-point function by 

(Tr X a X a (p J )x=i = -rj- 2) 2 {J - l)h? (Tr W(<p) Tr ^ J+2 ) conn . (2.14) 

The relation f)2.14p may be derived as follows^: The definition (12 .11) implies that the J\f = 1 
correlators are related to the M = 2 correlators by 



(0) M=1 = (0) M=2 + (0 [ a m A Tr W{T)) N=2 + ■ ■ 

Jc 2 



(2.15) 



Hence we have 

(/ «( ,2) A Tr J rJ+2 )x=i = hz 1 ^ 2 ^! ip J+2 } + h 2 z 1 z 2 z 1 z 2 (Tr (p J+2 TrW(<p)} + • • • , (2.16) 
Jc 2 

V(l) x=l = h 2 z 1 z 2 z 1 z 2 ((l) + hz 1 z 2 (TrW( V )} + --^ . (2.17) 

Note that in the expansion the first term vanishes after setting z 1 ,^ 2 —>■ 0, while the higher- 
order terms vanish because of z 1 ,^ 2 — > 0. Summing up all the terms in (12.141) into the 
exponential function, we find that the glueball loop operator is given as 

(Tr \ a \ a e wp )tf=i = — ^(Tr W((p) Tr ^) conn . (2.18) 

The more conventional quantity is the glueball resolvent, which is the Laplace transforma- 
tion of the glueball loop operator 

R(z) := (Tr X a X a — - — ) = [ du e" 2n (Tr X a X a e u *) . (2.19) 
z-<P Jo 

The explicit form of the glueball resolvent is known to be 



R( z ) = -W\z) + ^W'(z) 2 - Af(z) , (2.20) 

from the relation to the matrix model [4] or the generalized Konishi anomaly [11]. For the 
superpotential W(z) of degree n + 1, we need to include a polynomial f(z) of degree n — 1. 
Since we are considering the maximally confining phase, we have to choose the coefficients 
of f(z) so that ^W'(z) 2 — 4f(z) is factorized into a product of a polynomial and the square 
root of a quadratic function. Therefore, the resolvent R(z) in the maximally confining phase 
is characterized by the following conditions: 

• The resolvent R(z) is a linear combination of 1 and the square root of a quadratic 
function with the coefficients being polynomials of z. 

• The behavior of z —>■ oo is (Tr X a X a ) ■ z' 1 + 0(z~ 2 ). 

• The coefficient of 1 is —W'(z). 

We would like to reproduce all these features from the following instanton calculus. 



3 Two proofs of (|2.14j) have been presented recently in [HH]. 
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3 Scalar two-point functions 



In this section we compute the two-point correlation function of the loop operator: 



(Tre^Tre^) = ^ ]T 



oo r -i k 



k=0 



h 2 



Ch Y (a)Ch Y ([3) 



2 ' 



where 



Ch y (a) = l + sh 2 (a)^e 



a(c(0)) 



(3.1) 



(3.2) 



DEY 



with a — th and (3 = uh. (We set a = for simplicity in this section, but we can easily 
reproduce the contribution by multiplying the final result by the factors e ta and e ua .) Taking 
care of the constant term separately, we find 



(Tr e t(p Tr e uip ) = W + (Tr e tlp ) + (Tr e uip ) - 1 , 



where 



i 00 r i k 



fc=0 L 



|y|=k 



(3.3) 



(3.4) 



3.1 Operator formalism 

It is well-known that there is a correspondence between the Young diagrams and the fermion 
Fock states with neutral charge. By making use of the correspondence, we can compute the 
summations over the set of Young diagrams in operator formalism. Let us introduce a pair of 
charged (NS) free fermions 

${z) = £ , P{z) = £ * , (3.5) 

with the anti-commutation relation 

{A,^* s } = $r+s,o , r,seZ+^. (3.6) 

We define the Fock vacuum |0) by 

Vv|0> = V s *|0) = , r,s>0. (3.7) 

Recall that the Young diagram is specified by a partition A = (A*), where A; is the length of 
the i-th row. Then the corresponding state is given by 

OO 

l A > =II^-A i -iliO» , (3.8) 
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with -0* 1 0)) = 0, Vs. One can show that 

\x\=k 



J - J0> ^S.ira) |A) ' (3 ' 9) 



where J_i is the constant mode in the standard U(l) current J(z) :=: %p{z)ip*{z) :. From ( 13.91) 
we can confirm the following famous relation 

e n mo)- 2 = l , (3.io) 

\Y\=kO£Y 

which is used for computing the summation over the Young diagrams appearing in the partition 
function Z. We find (12.121) has a simple form: 

^ = exp(|) . (3.11) 

This is the most fundamental example of the computation in operator formalism of the sum- 
mation with the Plancherel measure [T2] . 

To compute the correlation functions, it is convenient to introduce the operator 



S^(z) :=: ^z)r(e^z) := ]T^— 1 , (3.12) 



which is a "point splitting" deformation of J(z) = E^(z). The modes of £( a '(z) 

4«> := e<r ~ %) ■ ^-r€ ■ , (3.13) 

rez+i 

satisfy the commutation relation^, 

[4 a) ,e£>] = MnP - ma) e$£> + 6n +m /^^® • (3.14) 
We will use the following important relation in the computation of correlation functions: 

oo 

£^ I A) = (e a ( x '- i+ ^) - e Q H + ^)) | A) = sh(«) J^exp (ac(n)) |A) . (3.15) 
i=l aex 

The second term comes from£,S a) ||0)) = -(sh(a))- 1 ||0)), which can be calculated directly from 
the definition of ||0)) or from the consistency £q |0) = 0. 



6 The infinite-dimensional Lie algebra with the commutation relation (|3. 14|) appeared first in |15j . where it 
was called area-preserving torus diffeomorphism algebra. The representation theory was initiated in [16] and 
the algebra was identified with W\+oo algebra. Recently the same algebra is used to reveal a connection of the 
counting of plane partitions and the Toda hierarchy [17] , 
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3.2 Computation 



After the preparation of the operator formalism, let us proceed to the explicit calculation of 
W defined by (13.41) . First of all, note that the formulas in the operator formalism imply 

(OlJ^^f Alo) = (fc!) 2 sh(a)sh(/3) Y, Dey n tmS • ( 3 - 16 ) 

Thus the problem of evaluation of W reduces to the computation of (0| J^SqSq J^jJO). With 
the help of the commutation relation (13.141) with J\ = S[°\ we obtain 

(q\ jk £ (<*) g (P) jk , v _ y- \p (fc!) 2 sh*(oQsh m (/3) , , (a) tf) l+m , . , , 

Z=0 m=0 v / V / 

where we have used 

jk-l-mjk^ = — ^JiriO) • (3.18) 

Bringing J_i's to the most left in (0| • ■ ■ |0) by iteratively using the commutation relation 

[S^elP, J_d = £[ a) [&W%-x + <W) + (sh(a)^g + 5 ltl )£<g) , (3.19) 

we are left with a constant term and terms such as (0|£„ |0). The coefficients of these 
terms can be calculated as follows. We first place the operator S^Sm ^ at the point (I, m) in 
the two-dimensional plane. (See Figure 1.) Using the commutation relation (13.191) once means 
that we bring the operator downwards or leftwards by one unit in the plane. After (I + m)- 
times of the iterative moves, we finally arrive at the integer point (n, —n). The coefficient of 
the (0|£n£^n |0) term (or the constant term) can be interpreted as the combinatorial factor 
of moving from (l,m) to (n, —n) (or (0,0) respectively) by these iterative moves. In this way 
we find 

(0jg/ a) gff Ji+ m lQ) [sh(a)} l - n [sh((3)} m+n sh n(a + (3) [shCa)]'- 1 ^)]"- 1 



{l + m)\ (I - n)\(m + n)\ sh (a + (3) l\m\ 



(3.20) 



Note that the second term, which comes from the constant term in the commutation rela- 
tion (13.191) . contributes only when I ^ and m ^ 0. Plugging back to (13 .4p . we call each 
contribution from the above two terms as Wi and W2, respectively. 

Let us concentrate on W2 first. 



" ' k ~' r V ' [sh(a)] 2l [sh(P)] 2m 



fc=0 1=0 m=0 



Q_ 

h 2 



(l\) 2 (m\) 2 (k-l-m)\ 



©Z>0@m>0 • (3-21) 
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Figure 1: Combinatorics of the coefficients. 



We can exchange the order of summation by bringing the fc-summation into the most right 

oo k k—l oo oo k—l oo oo oo 

£££•••=£££•=££ £ ••■ 

k=0 1=0 m=0 1=0 k=l m=0 1=0 m=0 k=l+m 

After the exchange of summation, we find W2 is computed to be 

W 2 = (/o(A)-l)(/ (5)-l) , (3.23) 
with A = 2y/qsh(a)/h, B = 2 v /gsh( ( 5)/ft and I n (x) is the n-th modified Bessel function. 
Now let us turn to W±: 



oo k k—l I 



Q_ 

h 2 



^. = i££££ 

fc=0 1=0 m=0 n=0 

Here the exchange of summation goes as 

00 k k—l I oo oo 00 I 



sh(a)] 2 *- n+1 [sh(/?)] 



2m+n+l 



shn(a + B) 



/!(/ — n)\m\(m + n)\(k — I — m)\ sh (a + 3) 



(3.24) 



00 / 00 00 



00 00 00 00 



££££■■■ = ££ £ £■■ = £££ £ -=£££ £ 



k=0 1=0 m=0 n=0 



1=0 m=0 k=l+m n=0 



1=0 n=0 m=0 k=l+r 



n=0 l=n m=0 k=l+m 



(3.25) 



where we have used ( 13 . 22 j) in the first equality and then brought the n-summation to the most 
left. Performing the /c-summation, we have 



sh (a + 3) . 

n=0 K y > x l=n 



Q_ 

h 2 



[ah(a)] 



21-n+l 



/!(/-n)! 



£ 



m=0 



/? 2 



[sh(/3)] 



2m+n+l 



which implies 



Wi = sh(ar)sh(/3)£; I n {A)I n {B) 



n=0 



m\(m + n)! 

sh n(a + /9) 
sh(a + 3) 



(3.26) 



(3.27) 
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Plugging back into (13. 4p all our results including the one-point function (Tre^) = Io(A), 
we finally find that the two-point function is given by 

(Tre^Tre^) = I (A)I (B) + sh(a)sh(/3) I n (A)I n (B) ah ^ a + ® . (3.28) 

n=0 S \ Q + IJ) 

Note that the first term can be interpreted as the disconnected contribution. Hence, 
(Tre^Tre<^) conn := (Tre^Tre uv ) - (Tr e ttp ) (Tr e utp ) 

= sh(a)sh(P)J2UA)In(B) 8 ^^^ . (3.29) 



ra=0 



4 Glueball one-point functions 



Having finished our computation of the scalar two-point function in the previous section, let 
us turn to the glueball one-point function. As (I2.18P implies, all the information we need is 
encoded in the scalar two-point function (I3.29P in the limit fi — > 0: 

1 oo 

lim -<Tre^Tre^) conn = e ta e ua tu ^ nl n (2y/qt)l n (2y/qu) , (4.1) 



n=0 



where we have restored the classical value a = (Trip) by simply multiplying the factors e ta 
and e ua . Using this result, we find only the exponential term in ( 12. 18ft gives a non-trivial 
contribution: 

(Tr \ a \ a e ULp ) = lim--^(TrW%)Tre^ conn . (4.2) 

Then, all we have to do is to pick up necessary terms out of (14. ip and perform the Laplace 
transformation to obtain the glueball resolvent. 

4.1 An example 

Before proceeding to the general superpotential, let us consider a simple example 

W(z) = 4- , (4-3) 



which allows an explicit computation and is instructive for general cases. Since the connected 
two-point function is given by 

lim — (Tr e tip Tr e ulp ) coaa = e ta e ua tu 

x (l x (2y/qt)I x (2jqu) + 2I 2 (2^qt)I 2 (2^qu) + 3I 3 (2jqt)I 3 (2jqu) + •■■), (4.4) 
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with the modified Bessel functions being 
we find 



0!1! 



1!2! 



>/2f 
0!2! 



(4.5) 



+ 



1! 



0!2! 



2ul 2 (2y/qu)e v 



+ 



(4.6) 



by picking up the t 4 terms. Therefore the one-point function of the glueball loop operator is 
given as 



(Tr\ a \ a e u *) = -2^ 



■a 3! 



9 r 2 31 



0!3! 
2! 2 



\{ 2 ^qu)e ua + ^--I x {2^qu)^ 



0!3!m 



l!2!u 



0!2!m 



h{2jq U y 



3! 



1! 1 



(4.7) 



which is transformed into the glueball resolvent through the Laplace transformation. The 
Laplace transformation can be performed as 



/ due-^ u -I n (2^u) = Z T 
Jo u 



(4- 



with Z (and Z which will appear later) defined by 
Z 



z — a 



\J{z - a) 2 - 4q — _ z - a + ^J(z - a) 2 - Aq 
, Z — 



2^q - 2^q 

Here we have used the Laplace transformation formula for the modified Bessel function 



(4.9) 



dte~ st I n (ujt) 



(s — ys 2 - lC 



2\n 



2,.,n 



to V s 2 — uj^u: 

and the recursive relation of the modified Bessel function 

2n 



x 



-I n (x) = I n -i{x) - I n+1 (x) 



(4.10) 



(4.11) 



Using (14.81) we easily find 
(TrA Q A Q -^) = -2 v ^ 3 ] 



z — ip 



0!3! 



0!3!' 



3! 

1!2! 



9 r 2 31 



2! 



2! 

0!2? 



2Jqa 2 — 
v 2!1! 



0!1! 



(4.12) 



It is remarkable that the Laplace transform of the modified Bessel function gives a linear 
combination of 1 and a/ \z — a) 2 — Aq with the coefficient being the polynomials of z, which 
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is required for the resolvent in the maximally confining phase. Furthermore, it behaves as 
0(z~ l ) in the limit z — > oo, as can be seen from 



Z = Z 1 



(4.13) 



Hence, in the final step of matching our computation to the expected result, all we have to 
do is to prove that the coefficient polynomial of 1 is — z 3 . This can be done without explicit 
calculation. Let us first rewrite our final result (I4.12p as 



V{z) + H{z)y/(z - a) 2 - Aq 



-2^f 



3! 

0!3! 



A Z 3 , A z 

0!3! 1!2! 



IJq a—— 
v 1!2! 



0!2! 



3! 



0!1! 



(4.14) 



Then, our task is to prove V(z) = z 3 . To pick up V(z) let us consider the "conjugate" of 

(OH) 



V{z) - H{z)\J{z- a) 2 -4g 



3! 

0!3! 



0!3! 



1!2! 



9 3! 



2 V 



0!2! 



3! 

2!l! 



11 r Z 



0!1! 



and add up with original (14.141) . Then we can sum up the right hand side into 



-2V{z) = -2^ — [Z + Zf - 2^fa 



1!2! 



2Vqa 2 ^- t \Z 
v 2!1! L 



(4.15) 



Z] , (4.16) 



because of (14.131) . We can further rewrite the result into 



-2V{z) = -2^f 



— a 
Z + Z+ — 



+ 2F 



with Fa defined as 



3! 2! a 3! 0! 

+ 



Using 



we find (14.171) reduces to 



l!2!l!l!yg 3!0! 0!0! 



— a z 
Z + Z + — = — 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



V(z) = z 3 -F 4 . 

To finish our proof of V(z) = z 3 we need to choose a = so that Fa_ vanishes. 

Note that since V(z) is the first derivative of the potential (in the matrix model terminol- 
ogy), we can interpret the constant part F 4 as the force felt by the cut. (See Figure 2.) Here 
the main contribution of _F 4 is _F 4 ~ a 3 for large a and the subleading term depends on the 
width of the cut y/q. 
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F 



Figure 2: Force felt by the cut. 
4.2 Monomial superpotentials 

The above argument can be easily generalized to any monomial superpotential. To compute 
the glueball one-point function for a superpotential 

7 k+l 

W(z) = ^ , (4.21) 
we have to pick up the t fc+1 -terms in the right hand side of (14.1 1) : 

fc-1 /, n [(*-I-l)/2] , _ , _ _ 

(Tr X a X a e u ' p ) = -2e ua ^ ^ k ~ l a l I ) £ I * 1 h-i-2 m (2^qu) . (4.22) 

1=0 ^ ' m=0 ^ 7 

Using the Laplace transformation formula (14.81) . we find the resolvent is given by 

fc-1 /r.\ Kfc-'- 1 )/ 2 ] 

7/ '-'- 2 '" (4.23) 



Z=0 v 7 m=0 v 7 



Note again that the result takes the form of —V(z) + H(z)^Jz 2 — 4q. To see the result is the 
expected one, we have to show V(z) = z k . This can be done by adding the "conjugate". In 
this way, we will complete "half" of the binomial expansion in (14.231) into a full one except 
constant terms. Finally we find V(z) = z k — Fk+i with the force Fk+i given by 

[fc/ 2 l 

If the potential is even, the force is odd and we can always make F k+ i = by setting a = 0. 
This means that in the even potential we can stabilize the cut at the center. However, this 
does not work for the odd monomial potentials. 

4.3 General superpotentials 

If the superpotential is W2(z) = z 2 /2, the glueball resolvent is given as 

1! 



(TrA a A Q — - — ) = -2y/q-^' 



z-ip' v ^0!l! 



0!1!' 



(4.25) 
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and the force felt by the cut is given as 



2 Vg Vl!0!0!0!^g 



while if the potential is W%(z) = z 3 /3, the glueball resolvent is 



(Tr A a A a ^— ^— ) = — 2 v /g 2 2 



z — (f 



0!2! 



0!2!' 



0!!!' 



and the force is 



2 i 2! 2! 2! 0! 
3 ~ ^ ' 0!2!1!1! + 2!0!0!0! 



Hence for the potential of the linear combination 
we have to impose the force balance condition: 



r / 1! 0! a \ ^ 2 /2! 2! 
mF 2 + gF 3 = mV?^^— j + ^ 



2! 2! 2! 0! 

+ 



2!0!0!0! 



. 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



As is clear from the above example the linearity holds for any superpotential. Therefore, 
as long as we are careful with the force balance condition, we can reproduce the results of all 
the superpotentials in the maximally confining phase. 



5 Higher genus correction 

So far we have computed the classical limit of the glueball one-point function by instanton 
calculus and found a complete agreement with the matrix model result. In the computation 
we have to introduce an equivariant parameter H for the toric action on C 2 , which can be 
related to the (local) SO (A) rotations of the ^-background [2]. 

Note that the /l-expansion of the scalar one-point and two-point functions computed in 
this paper agrees to the genus expansion of the Gromov-Witten theory of P 1 developed in 
[TBI [14] . The correspondence is described as follows: Okounkov and Pandharipande computed 
all genus correlation functions of the Kahler class uj and its descendents t p (u) (the so-called 
stable sector) with respect to the two ramification data at the north and the south poles of 
P 1 , which are labeled by the partitions. To obtain the gauge theory correlation functions we 
identify the operator Tr^ 2 - 7 as the cohomology class r p {uS). Then, the fc-instanton sector is 
recovered by taking both the ramification data to be (l fc ). 
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On the other hand, the same glueball one-point function was evaluated in pf] by relating 
the holomorphic quantities in supersymmetric gauge theory to amplitudes of topological string 
theory [10] and computing with the matrix model using open/closed string duality. Here 
the loop expansion parameter of the matrix model was interpreted as the genus expansion 
parameter of topological string theory, which, on the gauge theory side, was identified with 
a constant graviphoton background. A comparison between instanton calculus and matrix 
model at higher genus through topological string amplitudes was made in [TSJ [T9] . 

Although it was pointed out [20J that the f2-background and the graviphoton background 
are different^, here we would like to compute the first order correction to the one-point function 
on both backgrounds and compare with each other, because these two backgrounds look 
similar and share the interpretation of genus expansion. We shall start with computing the 
^-correction of the glueball one-point function from the instanton calculus in subsection 5.1 
and proceed to recapitulating the matrix model computation in subsection 5.2. Comparing 
the result of the glueball one-point function from the instanton calculus (15.41) with the result 
of the one-point function in the matrix model (15. lip or (I5.14p . we find a discrepancy. Though 
we cannot find any way to relate the two results, we still expect there to be a relation between 
these two genus expansions, for example, by change of variables. 

5.1 /i-corrections of the glueball one-point function 

Let us study the ^-corrections of the glueball one-point function. We shall choose the simplest 
Gaussian superpotential here. 

Since the connected two-point function is given as 




(5.1) 



we find the ^-corrections of the glueball one-point function is 




(uh) 2 



24 



) 





(5.2) 



where we have used identities of the modified Bessel function 



2n 



I n (x) , + I n +i(x) = 2—I n (x) . 



(5.3) 



x 



Therefore the glueball resolvent is 




[{z - a) 2 - 4g] 5 /2 



2(z-a) 2 + Aq 



(5.4) 



7 The relation between these two backgrounds was also discussed in [21] . 
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Note that the following formula holds for the Laplace transformation 

C{^f(t)} = (-l) n F^(a) , (5.5) 
if we denote the Laplace transformation as £{/(£)} = F(s). 

5.2 Matrix model computation 

The computation on the matrix model side can be performed with several methods. The first 
method is due to loop equation: The resolvent in the matrix model is determined by the loop 
equation. 

We shall consider the matrix model with potential W(M): 



Z = J VMexp^-NTiW(M)^ 



(5.6) 



The loop equation for the resolvent R(z) = (N 1 Tr (z — M) 1 ), 
(See (13.52) in [22\.) gives 

at one loop. For the gaussian model W(y) = fiy, the one-loop contribution is only non- 
vanishing for (Tr (p n ) with n > 4. Hence, the resolvent for z — > oo should behave as R\{z) ~ 
1/z 5 . This means there is no pole at y = oo and we can perform the y integration easily by 
picking up the pole at y = z: 

(2R (z) - W'{z))R x {z) = -JL^Rv{z) • (5.9) 

The right-hand-side is a two-point loop correlator and given by [23] 

Ro(z) = R (z, z) = . - 1 A . . (5.10) 



SW(z) uw uv ' ' fi(z 2 -4/fi) 2 ' 
Therefore we find 

Rl{z) = tf{z* - 4//^) 5/2 ' (5 ' H) 

For the case of the Gaussian matrix model, an exact manipulation is possible. The exact 
solution for Gaussian matrix model 

<^Tre^> = iy PM^Tre^ exp(-^TrM 2 ) , (5.12) 
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is given by [24] 



(_Tre uM ) 

\ AT I 



u 



h(2u/y/Jl) + 




(5.13) 



up to genus one. After the Laplace transformation we find 




z-M 



1 



N 2 ^ 2 (z 2 - A/fi) 5 / 2 



1 1 



(5.14) 
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